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Abstract. Totally asymmetric simple exclusion processes on lattices with junctions, 
where particles interact with hard-core exclusion and move on parallel lattice branches 
that at the junction combine into a single lattice segment, are investigated. A simple 
approximate theory, that treats the correlations around the junction position in a 
mean- field fashion, is developed in order to calculate stationary particle currents, 
density profiles and a phase diagram. It is shown that there are three possible 
stationary phases depending on the state of each of the lattice branch. At first- 
order phase boundaries, where the density correlations are important, a modified 
phcnomenological domain-wall theory, that accounts for correlations, is introduced. 
Extensive Monte Carlo computer simulations are performed to investigate the system, 
and it is found that they are in excellent agreement with theoretical predictions. 
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1. Introduction 

Asymmetric simple exclusion processes (ASEP) have been introduced originally in 1968 
as a theoretical model for a description of kinetics of biopolymerization pQ . Although in 
last years the area of application of ASEP has been significantly broadened j2j OH Ej , and 
it now includes a road traffic flow analysis IE] , polymer dynamics in dense medium 
[7j, and many other problems, the major application of ASEPs remains the modeling of 
various biophysical transport phenomena. In particular, simple exclusion processes have 
been used successfully to describe protein synthesis jHUH], mRNA translation phenomena 
[Til] , gel electrophoresis ^JJ , a motion of motor proteins along the cytoskeletal filaments 
[T2] and the depolymerization of microtubules by special enzymes [T3] . 

ASEP are discrete non-equilibrium models that describe stochastic dynamics of 
multi-particle transport along one-dimensional lattices. The lattices are finite and 
generally consist of L ^> 1 sites. Each site can be occupied by a single particle or 
empty, and the particles interact only through the hard-core exclusion potential. The 
dynamics of ASEP is asymmetric, i.e., the particles can hop to the left or to the right 
but with different probabilities. In the simplest totally asymmetric simple exclusion 
process (TASEP) the particles move only in one direction. In this case the rules for the 
motion are the following. The particle at site 1 < i < L can move one step forward if 
the site i + 1 is empty. The particle can enter the lattice with the rate a if the first 
site is available, and it can leave the system from the last site i = L with the rate j3. 
In the stationary-state limit of TASEP, the system can be found in one of three phases 
depending on entrance, exit or bulk processes dominate the overall dynamics. 

The unusual dynamic properties and phase behavior of ASEP and a wide range of 
applications in chemistry, physics and biology stimulated many theoretical studies of 
asymmetric exclusion processes [2 El El- There are several exact results for the steady- 
state properties of ASEP for different update rules [31 E], although most theoretical 
investigations utilize the approximate methods along with Monte Carlo computer 
simulations [U EH EH EH]. The coupling of several exclusion processes has been considered 
in the study of parallel-chain ASEP [IB*] ITT] ITS] The combining of non-equilibrium 
exclusion processes with equilibrium particle association/dissociation phenomena led to 
unusual phenomena of localizations of density shocks PHI EH1 EU 1221 E3L However, 
despite the differences in the specific dynamic rules and stationary properties of specific 
ASEP, the microscopic origin of unusual dynamic properties and phase behavior can be 
well understood with the help of a phenomenological domain wall theory [2U E] . 

The majority of the investigated asymmetric exclusion processes deal with the 
particle movement along the one-channel lattices. Although the one-channel approach 
describes many situations in the biophysical processes, the more realistic description of 
the cellular transport requires an extension of original ASEP to include the possibility 
of a transport on the lattices with a more complex geometry. For example, consider 
the motor proteins kinesins that move vesicles and organelles along the microtubules 
and play important role in cellular transport [221 E] • Microtubules are made of parallel 
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linear polymers, called protofilaments, that are arranged circumferentially. It is known 
that kinesins walk only on single protofilaments. However, experiments [2H| indicate 
that the number of protofilaments may vary, at least for in vitro conditions, and it 
indicates the existence of junctions and other lattice defects. Such defects might lead to 
motor proteins crowding phenomena that are responsible for many human diseases j27j . 
These observations suggest an importance of investigation of the asymmetric exclusion 
processes on the lattices with junctions as a model for these complex biological transport 
phenomena. 

Recently Brankov et al. [2H] have investigated TASEP on chains with a double- 
chain section in the middle by using an approximate theory and computer simulations. 
This corresponds to having two consecutive junctions on the lattice. Several stationary 
phases, the existence of which has not been expected, are found. In addition, the density 
profiles at phase boundaries and strong correlations between different lattice branches 
have been observed, but they could not be explained because the theory neglected 
the correlations. In this paper, we present a theoretical investigation of TASEP on 
the lattice with one junction. We develop a simple approximate theory that allows 
to calculate particle currents, density profiles and a phase diagram at large times. At 
first-order phase boundaries, we develop an extension of the domain wall theory [21] 
that accounts for correlations. This approach is then applied for explicit calculations of 
density profiles. The theoretical predictions are compared with extensive Monte Carlo 
computer simulations. 

The paper is organized as follows. In section 2 the model is introduced and 
theoretical calculations of stationary properties are presented. In section 3 theoretical 
predictions are compared and discussed with the results of Monte Carlo computer 
simulations. The final summary and conclusions are given in section 4. 

2. Theoretical Description 

2.1. Model 

We consider identical particles that move along the lattice with the junction positioned 
in the middle of the system as shown in Fig. 1. The system is out of equilibrium, and 
it has three equal-size branches, each containing L sites. The particles can enter chain 
I or chain II with the rate a if the first site at the corresponding branch is available. 
Two chains merge together at site L + 1 and form chain III. The particle can leave the 
system with the rate (3: see Fig. 1. Inside of the lattice chains the particle can only 
move one step forward if the neighboring site is empty. 

Without chain I or chain II, the junction disappears and the system reduces to a 
totally asymmetric simple exclusion process on the one-channel lattice, for which the 
full description of stationary properties, such as phase diagram, particle currents and 
density profiles, are known El H]- This simplest model of TASEP on the lattice 
without junctions has three stationary phases. When the entrance into the system is a 
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Figure 1. Schematic picture for the model for TASEP on the lattice with a 
junction. Particles can enter chain I or chain II with the rate a. At the site 
L + 1 two lattice branches coalesce into chain III, from which particles can exit 
with the rate (3. Arrows indicate the allowed transitions, while crossed arrows 
correspond to the prohibited moves. 



rate-limiting process, for a < 1/2 and a < (3, the system is found in a low-density (LD) 
phase with the current and bulk density given by 

Jld — a(l - a), Pbuik,LD = a. (1) 

If exit controls the dynamics of the system, for (3 < 1/2 and (3 < a, the stable stationary 
state for the system is a high-density (HD) phase with the following current and bulk 
density 

JhD = P{1~P), Pbulk,HD = 1 - (3- (2) 

Finally, for large entrance and exit rates (a > 1/2 and f3 > 1/2), when the dynamics is 
determined by bulk processes, the system is in a maximal- current (MC) phase with 

JmC = 1/4, Pbulk,AW — 1/2- (3) 

Full density profiles for systems of any size can be calculated explicitly [2|. 

The microscopic origin of complex phase behavior and unusual dynamic properties 
of asymmetric exclusion processes can be explained via a phenomenological domain-wall 
(DW) theory [23]. According to this approach, the domain wall is the boundary region 
between two possible stationary phases, and it moves through the system as a random 
walker with a speed determined by the currents and densities in two phases: 

v DW = u + -u_ = — — —, (4) 
P+-P- 

where "+" ("— ") corresponds to the phase to the right (left) of the domain wall, and u + 
and m_ give the domain wall rates to hop to the right or left. For vdw > (u + > uJ) 
the domain wall moves to the right and the "negative" phase becomes a stationary state 
of the system, while for v^w < (u + < uJ) the domain wall travels to the left and 
the "positive" phase wins over. On the phase boundaries the domain wall has equal 
probability to go forward or backward, i.e., u + = w_ and vdw = 0. As a result, the 
density profiles are linear. This is due to the fact that the domain wall can be found 
with equal probability at any position in the system. 
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2.2. Theoretical calculations for stationary phases 

The overall state of the system is specified by the nature of phases that might exist in 
each of the lattice branches. Since only three stationary phases can be found in each 
lattice chain (HD, LD or MC), the total number of possible stationary phases in the 
system with the junction is equal to 3 3 = 27. However, because of a symmetry, chains 
I and II should have identical phases (see Fig. 1), and the overall state of the system 
is determined by phases in the chain I (or II) and the chain III. Then the number of 
possible stationary states reduces to only 3 2 = 9. 

The overall stationary current passing through the system can be written as 

Joverall = J III = -h + -hi = 2Jl- (5) 

It suggests that chain I and II cannot have the the maximal-current phase with J = 1/4, 
because the maximal possible current through the system is just equal to 1/4. Thus, 
there are only 6 possible stationary phases: (LD,LD), (LD,HD), (LD,MC), (HD,LD), 
(HD,HD) and (HD,MC), where in the expression (A,B) A describes the phase in the 
chains I and II, while B corresponds to the phase in the chain III. 

The junction introduces an inhomogeneity in the system and it makes impossible to 
solve the large-time dynamics exactly ^5EH1- However, TASEP on the lattice with the 
junction can be mapped into 3 coupled homogeneous asymmetric exclusion processes 
as shown in Fig. 2, for which the approximate description can be developed [T%} I29j. 
In order to obtain the dynamic properties of the system explicitly, in the simplest 
approximation, we assume that there are no correlations in the occupation of the sites 
before and after the junction, i.e., 

Joverall = ^junction = 2 < T L (1 - T L+1 ) >^ 2 < T L > (1- < T L+1 >), (6) 

where < t l >= p L is the probability to occupy the site L on the chain I or II, and 
< tl+i >= Pl+i is the average particle density at the site L + 1 of the chain III. The 
effective rates a e ff and j3 e ff (see Fig. 2) can be expressed in terms of the particle 
densities at the sites near the junction, 

OCeff = 2pL, Peff = 1 ~ PL+1- (7) 

Now we can investigate the existence of different stationary phases. Consider first 
(LD,LD) phase, which can be specified by the following conditions, 

a < 1/2, a < p e f f ; a eff < 1/2, a eff < p. (8) 

The stationary currents and bulk densities are given by 

Jj = a(l-a), pi,butk = at, Jni = « e //(l - a e //), Pi n, bulk = a e//-(9) 

Using the expression (0) for currents, the effective rate a e ff can be expressed in terms 
of the entrance rate a, 



— 8a (1 — a) 

OL e ff = * . (10) 
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Figure 2. The totally asymmetric simple exclusion process on the lattice with 
the junction can be mapped into three homogeneous TASEP coupled at the 
ends. The effective rate /3 e // describes the process of exiting from the chains I 
or II, while the effective rate ot e ff corresponds to the entrance into the chain 
III. 



Thus a e ff always satisfy the conditions (JBJ), however, this equation yields physically 
reasonable values of the effective entrance rate only when the term in the square root 
is positive, i.e., 

1 \/2 

a<-- — « 0.146. (11) 

2 4 v ' 

Because Pl+i = a e // 111 (LD,LD) phase, the effective exit rate, as follows from Eq. (J7J), 
is given by (3 e jf = 1 — a e ff- Then the condition a < (3 e ff is true for all values of 
parameters. Thus the system is in (LD,LD) phase when 



1 - ,/l - 8a(l - a) 1 v/2 

(3 > * and a < - - — . (12) 

1 2 2 4 v ' 

For the (LD,HD) phase the conditions of existence can be written as: 

a < 1/2, a < (3 eff ; (3 < 1/2, (3 < a eff ; (13) 

while the stationary currents and bulk densities are 

J/ = a(l-a), pi,buik = a, J in = (3(1 - (3), pm^ik = 1 - (3- (14) 

Because the particle current is stationary, Eq. (j3j) implies that 



1 - ,/l - 8a(l - a) 1 J2 

(3 = ± -, and a < - - — . (15) 

2 2 4 v ; 

These expressions describe the parameter's space for (LD,HD) phase. 

Similar analysis can be performed for (LD,MC) phase. The allowed parameters for 

this phase are specified by 

a < 1/2, a < (3 eff - a eff > 1/2, (3 > 1/2. (16) 
The currents and bulk densities have the following values: 

J I = a(l-a), pi,buik = a, Jin = 1/4, pm,buik = 1/2. (17) 
Then the currents in the chain I and II are Jj = Jn = 1/8, which implies that 
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This equation along with the condition (3 > 1/2 fully determines the region of existence 
of (LD,MC) phase. 

The situation is very different for (HD,LD) phase for which 

/3 eff < 1/2, a > (3 eff ; a eff < 1/2, a eff < (3. (19) 

The corresponding equations for the currents and bulk densities are 

Jl — fteff(l—(3eff), Pl.bulk = l~Peff, Jill = a eff(l — a eff), PlII,bulk = 

It is known that in this phase pi, — 1 — Peff and pl+i = cteff- Then from Eq. (JZj) it can 
be shown that 

a e f f = 2(1 - p eff ), aadp e ff = l-a e ff. (21) 

However, these two equations have no real solutions together, and therefore (HD,LD) 
phase cannot exist for any value of the entrance rate a and the exit rate f3. 
The (HD,HD) phase is determined from the conditions 

Peff<l/2, Pe ff <a; /3 < 1/2, (3 < a eff] (22) 

The stationary properties of this phase are given by 

Jl — Peff(l—Peff), Pl,bulk = ^~Peff, J III = P{1— f3), PlII,bulk = 1~P- 

The stationary condition for the particle currents [see Eq.©] helps to determine the 
effective exit rate constant, 

^ WEWZg . (24) 

After combining this result with the set of phase existence requirements (J22j) we obtain 
the final conditions for (HD,HD) phase: 



P<-, if a > ^ ; (25) 

and 

P < 1/2, if a < 1/2. (26) 

The last possible phase in the system is (HD,MC) phase, which is specified by the 
following conditions: 

Peff<l/2, P eff <a; P>l/2, a eff >l/2. (27) 

The particle currents and bulk densities in this phase are given by 

Jl — Peff(l—Peff), Pl,bulk — 1—Peff, JlII = l/^, PlII,bulk = 1/2.(28) 

From Eq.© it can be easily shown that 

Peff = \ ~ ^- (29) 

Comparing this result with the conditions (|2*7j) . we finally derive 

a>i-^, and/3>^. (30) 
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These expressions describe the phase space for (HD,MC) phase. 

Thus the analysis that considers TASEP on the lattices with junctions as coupled 
asymmetric exclusions systems suggests the existence of 5 phases. The calculated phase 
diagram is shown in Fig. 3. It can be seen that two phases, namely, (LD,MC) and 
(LD,HD), correspond to stationary phase boundaries. This result is due to the neglect 
of correlations around the junction point. However, at these conditions, the domain 
wall, that separates two phases can move with equal probability to the right or to the 
left j21] • Then the density profiles will be different from the one predicted in the simplest 
theory, that treats the correlations around the junction in a mean-field manner. It can 
be expected that the density profile is linear in the chain I and the chain II segments for 
(LD,MC) phase, extrapolating between coexisting low-density and high-density phases, 
while the density profile in the chain III still corresponds to the maximal-current phase. 
For (LD,HD) phase the linear density profiles can be found in all lattice segments, 
although the slopes strongly depend on inter-segment density correlations. 
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Figure 3. Phase diagram for the totally asymmetric exclusion system on the 
lattice with the junction. Lines are calculated from theoretical predictions for 
phase boundaries, while symbols are from Monte Carlo computer simulations. 
Thick solid line corresponds to a non-equilibrium first-order phase transition 
between (LD,LD) and (HD,HD) phases. Thin solid line shows a continuous 
phase transition between (HD,HD) and (HD,MC) phases. Meanwhile, thick 
dashed line describes a mixed phase transition between (LD,LD) and (HD,MC) 
phases: a first-order transformation in the chain I and II, and the continuous 
phase transition in the chain III. The size of the symbols reflect the statistical 
error of computer simulations. 
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2.3. Phase boundaries 

Qualitative domain-wall arguments presented above indicate that particle density 
changes linearly for some phase coexistence lines. However, exact linear density profiles 
might differ significantly from the densities obtained within the simplest approximate 
theory that neglects the correlations near the junction. Here we utilize the domain-wall 
approach [21] to account for these correlations in order to derive the phase boundary 
density profiles explicitly. 

The line specified by a = \ — ^ for (3 > 1/2 (see Fig. 3) describes the mixed 
phase coexistence between (LD,LD) and (HD,MC). When crossing this line, the density 
jumps in the lattice segments I and II, while the change is continuous in the chain III. 
At the phase boundary the domain wall in the chain I or II moves randomly with equal 
forward and backward rates between the low-density (pld = \ ~ 4 ) and high-density 
(Phd = 1 — Peff = \ + 4 ) regions. The domain wall picture cannot be used in the 
lattice chain III when the maximal-current phase appears [21|, and the inter-segment 
correlations are not significant. Then the resulting density profiles in the left lattice 
segments are expected to be linear, connecting the low-density and high-density values. 

To determine the density profiles at the phase boundary between (LD,LD) and 
(HD,HD) phases is a much more complicated problem because of the strong correlations 
around the junction. This phase coexistence line is given by the following conditions: 



as shown in Fig. 3. The domain wall separates the coexisting stationary phases and it 
can be found in any lattice segment. Let us define a position of the domain wall in the 
system via relative coordinate x, 



where i is the site index and L is the length of one lattice segment. Therefore, the case 
of < x < 1 describes the chain I and II, while 1 < x < 2 corresponds to the chain III. 

When the domain wall is in the chain I (and simultaneously in the chain II) it can 
move with the same rate ui to the left or to the right, while in the chain III it travels 
forward or backward with the rate unf. see Fig. 4. These rates can be determined by 
utilizing the expression (JH), 



a < - — 




(31) 



(32) 



Uk = 



P k +~P' 



for k = I, III 



(33) 



where 




1 — a 



P 



A P 



1-/3, 



(34) 



and 



Ji = a(l — a), Jin = 2Ji. (35) 
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As a result, we obtain the following expressions for the rates uj and um, 
a(l — a) 



Ui 



I -2a ' 
2a{\ - a) 



2a{\ - a) 



1 " 2/5 v /l-8a(l-a)' 



(36) 
(37) 




Figure 4. Schematic picture of the domain wall dynamics at the phase 
coexistence line between (LD,LD) and (HD,HD) phases. The domain wall 
in the chain I and II hops to the right or left with the rate ui, while in the 
chain III it moves with the rate um. 

To calculate the density profiles we introduce a probability Pj to find the domain 
wall at any position in the chain I or II, and Pju gives the probability that the domain 
wall is in the lattice chain III. These probabilities are obviously normalized, 

Pi + Pui = 1. (38) 

The probability that the domain wall occupies a specific site i is equal to Pi/ L or Pm/ L 
for i < L and i > L, respectively. Then at the junction 

ujPj/L = umPm/L. (39) 

This relation reflects the fact that the domain wall has equal probability to travel 
between different lattice segments. By combining last two equations we obtain: 

P, = ^-, P,„ = -^—. (40) 
ui + u IH ui + u IH 

These expressions have a simple physical explanation. The domain wall spends less 

time in the lattice segments where it fluctuates faster. All lattice chains have the same 

length L, and the domain wall in the lattice segment where it fluctuates faster is able to 

diffuse to the junction point quicker. As a result, it will jump to another lattice chain 



TASEP with Junction 



11 



more frequently. This is a reason for strong inter-segment density correlations at the 
first-order phase transitions. 

We can determine the probabilities of having the domain wall at any position less 
then a certain value of x. If the domain wall is in the channels I and II at the coordinate 
xdw, then this probability is given by 

Prob(x DW < x) = Pjx, < x < 1, (41) 
Similarly, for the domain wall in the lattice chain III, 

Prob(x DW < x) = Pi + P IH {x -l),l<x<2. (42) 
Then the density at any position can be calculated as 

p(x) = p k _Prob(x DW > x) + p k + Prob(x DW < x), k — I, IIP (43) 
Finally, combining Eqs. (j^ . (J5fij) . <|37 |) . (j^ . <fTTj> and we obtain 

2(1 - 2a) 2 

p(x) j = a H — , =x, < x < 1; (44) 

2(1 - 2a) + yT- 8a(l - a) 

and 



1-Jl- 8a(l - a) 2(1 - 2a) Jl - 8a(l - a) 
p(x)m = 



2 2(1 - 2a) + Jl - 8a(l - a 



l-8a(l-a) . 

+ 7 (45) 

2(1 -2a) + ,/l-8a(l-a) 



for 1 < x < 2. At the entrance sites we have, as expected, p(x = 0)i = a, while at the 

ie , 

2(1 -2a) 2 



last site p(x = 2)m = 1+ v /l 8 Q ( j_^ _ ^ _ ^ ^ ^ e junction, the densities are equal to 



2(1 - 2a) + Jl - 8a(l - a) 



p(x = l)i = a + 



l_.^l_8a(l-a) 2(l-2a)^/l-8a(l-a) 

p(x = l) m = s + = =. (46) 

2 2(1 - 2a) + Jl - 8a(l - a) 



It is important to note that the densities at junction are not equal to the values 1 — a 
for the chain I and II, and f3 = 1 V 1 a ) jjj^ respectively, as expected 

from the simplest approximate theory that assumes an independent coupling between 
the lattice segments. The domain-wall approach allows to take into consideration the 
correlations in the densities around the junction. 

3. Monte-Carlo simulations and discussions 

In order to check the validity of our approximate theory we performed extensive 
computer Monte Carlo simulations. Since the computer calculations for the cases of 
phase transitions are very time-consuming, especially for small values of entrance and 
exit rates, we utilized one of the continuous-time Monte Carlo algorithms, the so-called 
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BKL algorithm, first introduced by Bortz, Kalos and Lebowitz almost 30 years ago [30j. 
The main idea of the BKL algorithm is to create an event-driven update scheme, so 
that the rejected, "eventless" steps are skipped. The method is simple to implement, 
and with the latest improvements [3T| it is also very efficient and fast. 

In our simulations the number of effective steps per site was typically around 10 7 . 
At phase transitions we ran simulations much longer, and the number of effective steps 
per site was 10 8 -10 9 . For all simulations we neglect first 3% of Monte Carlo steps to 
account for the time that the system takes to achieve a stationary state. Our theoretical 
calculations assume infinite size lattice segments, however, in our simulations we used 
L = 100 and we checked that for larger sizes of lattice segments the results do not 
deviate from the ones presented here. 

A phase diagram obtained from Monte Carlo simulations is shown in Fig. 3. The 
boundaries between the stationary phases have been determined by considering the 
saturation in the currents and comparing qualitative changes in the density profiles. 
Specifically, the phase coexistence line between (LD,LD) and (HD,HD) phases is found 
when the density profiles become linear. The boundary between (LD,LD) and (HD,MC) 
phase is determined when the linear density profile is observed in the chains I and II, and 
the overall particle current saturates. Similarly, the current saturation method allows 
to specify the phase boundary between (HD,HD) and (HD,MC) phases. The overall 
error in the determination of phase boundaries is less than 5% ^H]- To illustrate our 
approach, the dependence of current on the entrance rate a (for fixed exit rate (3 = 1) 
is plotted in Fig. 5. It can be clearly seen that the current saturates at a ~ 0.15. 
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Figure 5. The particle current as a function of the entrance rate a. The exit 
rate is fixed, (3 = 1. The current saturates at a » 0.15 that corresponds to 
the phase transition between (LD, LD) and (HD, MC) phases. A solid line 
describes theoretical predictions, while symbols are the results of computer 
simulations. 

Monte Carlo simulations allowed also to calculate explicitly the particle densities. 
The resulting density profiles for different stationary phases are shown in Fig. 6. Since 
the stationary properties in the chains I and II are essentially the same, we investigated 
in detail only one of two equivalent lattice segments. The computer simulations for the 
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Figure 6. Density profiles for different bulk stationary phases: (a) (LD,LD) phase 
with (3—1 and three different entrance rates a = 0.05, 0.07 and 0.09; (b) (HD.HD) 
phase with a — 0.8 and three different exit rates (3 = 0.32, 0.38 and 0.44; (c) 
(HD,MC) phase with a = 1 and (3 — 0.58. Lines are our theoretical predictions 
for the full-length density profiles calculated from a, /3 e //, a e // and (3 by using the 
exact expressions derived in Ref. (0)- Symbols are obtained from the Monte Carlo 
computer simulations. 

densities in the bulk phases are well described by the theoretical predictions. Note, 
however, the deviations near the junction that grow as the system approaches the phase 
boundaries. 

The situation is very different for the phase coexistence line between (LD,LD) and 
(HD,HD) phase. Our theoretical approach predicts linear density profiles in all lattice 
segments, although with different slopes. The extensive computer simulations (see Fig. 
7) mainly confirm these suggestions, although the density profiles in the chains I and 
11 deviate slightly from the linearity. There are possible sources of these deviations: 1) 
it might be due to the errors in the determination of the exact position of this phase 
boundary; and/or 2) the cross-correlation between the particles in the chains I and II 
as was observed for the related system |2H] . 

The analysis of stationary properties of TASEP on the lattice with junction 
indicates an excellent agreement with the theoretical predictions: see Figs. 3, 5, 6 and 7. 
Although the simplest approximate theory neglects the correlations around the junction, 
it does not strongly affect the position of phase boundaries, stationary currents and bulk 
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Figure 7. Density profiles for phase coexistence line between (LD,LD) and 
(HD,HD) phases for different parameters: (a) a = 0.02, (3 = 0.0408; (b) 
a = 0.12, (3 = 0.3020. Solid lines are theoretical predictions from the modified 
domain- wall theory, namely, Eqs. (|44j) and (|45j) . Dashed lines correspond to 
the predictions of the simplest theory without correlations. 



density profiles. However, the effect of correlations is important at the first-order phase 
transitions between (LD,LD) and (HD,HD) phases. The extension of the domain- wall 
phenomenological approach, that argues that the domain wall fluctuates with different 
rates in the different lattice segments, is able to account for density correlations as 
compared with the results from the computer simulations (Fig. 7). 

4. Summary and conclusions 

The stationary properties of totally asymmetric simple exclusion processes on the lattices 
with junctions are investigated with the help of the simple approximate theory and by 
analyzing extensive computer Monte Carlo simulations. It is found that the phase 
diagram of the system consists of three stationary phases. This behavior is similar to 
asymmetric exclusion processes on the lattices without junctions, although the maximal- 
current phase cannot be sustained in the lattice branches before the junction due to the 
stationary current limitations. 

There are three different types of phase transitions in this system. On the 
coexistence line between (LD,LD) and (HD,HD) phase there are density jumps in all 
lattice segments that correspond to the first-order phase transitions. On the phase 
boundary between (HD,HD) and (HD,MC) the density profile changes continually in 
the lattice chain III. However, the phase transition between (LD,LD) and (HD,MC) has 
quite an unusual mixed character: the first-order transformation in the lattice chain I 
and II, and the continuous change in the lattice chain III. This phase behavior differs 
significantly from TASEP on the lattices without junctions. 

To analyze the dynamics of asymmetric exclusion processes on the lattice with the 
junction the approximate theoretical approach has been developed. According to this 
theory the system of particles moving on the lattice with the junction can be viewed 
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as three TASEP on the lattices without defects that coupled at the junction. As the 
simplest approximation, the correlations near the junction position are neglected. This 
method allows to calculate explicitly all stationary properties and the phase diagram. 
It is found that the theoretical predictions are in excellent agreement with Monte Carlo 
computer simulations for all phase regions except for the phase coexistence line between 
(LD,LD) and (HD,HD) phase. 

For the first-order phase boundary, where the inter-segment density correlations are 
important, the modified domain- wall approach is developed. We argue that the domain 
wall, that separates the low-density and high-density phases in each lattice segment, 
fluctuates with the different rates in the different lattice branches. It means that the 
domain wall does not spend the same time in all lattice segments, and this leads to 
the correlations observed in the system. The computer Monte Carlo simulations fully 
support the predictions from the modified domain-wall theory. It is suggested that this 
approach is general enough to be used successfully to account for correlations in other 
inhomogeneous asymmetric exclusion processes. For example, for TASEP on the lattices 
with double-chain sections [2E], i-e., the system with two junctions, can explicitly predict 
the density profiles and it can explain the observed density correlations. 

There are several extensions of the system that can be explored in a future. In the 
original model, the particle dynamics in the lattice segments before the junction have 
been identical. It will be interesting to investigate a system where two lattice branches 
before the junction are dynamically different. More complex dynamics is expected for 
a system where more than three lattice branches are joined together at the junction 
point. It is expected that our approximate theory without correlations and the modified 
domain-wall approach, supported by computer Monte Carlo simulations, will provide a 
reasonable way of analyzing these complex non-equilibrium systems. 
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